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Here we discuss the annealing behavior of an infinite-range ±J Ising spin glass in presence of 
a transverse field using a zero-temperature quantum Monte Carlo. Within the simulation scheme, 
we demonstrate that quantum annealing not only helps finding the ground state of a classical spin 
glass, but can also help simulating the ground state of a quantum spin glass, in particularly, when 
the transverse field is low, much more efficiently. 
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Quantum annealing (QA) 10 is a method of finding 
the ground state (minimum energy state) of a given clas- 
sical Hamiltonian employing external quantum fluctua- 
tions and subsequent adiabatic reduction of them. In 
QA, one is given with a complex classical Hamiltonian TL 
whose ground state is to be determined. In order to intro- 
duce the quantum fluctuations (necessary for the anneal- 
ing) into the system, one adds a quantum kinetic term 
TL'(t) to the classical Hamiltonian, such that TL'(t) and 
TL do not commute. Initially, one keeps TL'(t = 0) 3> TL 
so that the total Hamiltonian TCtot(t) — TL'(t) + TL is well 
approximated by the kinetic part only (TLtot(0) ~ W(0)). 
The system is initially prepared to be in the ground state 
of TL'(0) (one chooses TC'(0) to have a easily realizable 
ground state). Now since TLtot{0) ~ TL'(0), the overlap 
I (■)/>(£) l^i (i)} I between the lowest eigen- value state (we 
will call it instantaneous ground state) \<fii(t)) of the total 
Hamiltonian TLtot(t) and the instantaneous state \tp(t)) of 
the evolving system will be close to unity at t = (since 
V'(O)) is taken to be the ground state of TL'(0)). If one 
subsequently reduces TL'(t) slowly enough, then accord- 
ing to the adiabatic theorem of quantum mechanics, the 
overlap \(ip(t)\<fii(t))\ will always stay close to its initial 
value (i.e., unity). Hence at the end of such an evolution, 
when TC'(t) is reduced to zero at t = r, the system will 
be found in a state \tp{r)) with \{i(j(t)\4>i(t))\ « 1, where 
\4>x{r)) is the ground state of 7Y tot (r), which is nothing 
but the surviving classical part TL of the Hamiltonian. 
Thus at the end of an adiabatic annealing the system is 
found in the ground state of the classical Hamiltonian 
with a high probability. Based on this principle, algo- 
rithms can be framed to anneal complex physical systems 
like spin glasses as well as the objective functions of hard 
combinatorial optimization problems (like the traveling 
salesman problem or TSP) mapped to glass like Hamil- 
tonians, towards their ground (optimal) states. 



Here we study the annealing behavior of an infinite- 
range ±J Ising spin glass in a transverse field, using a 
zero-temperature quantum Monte Carlo. The model and 
the basic QA scheme for it are introduced in Sec|TJ In Sec. 
im we discuss at length the zero-temperature quantum 
Monte Carlo method used here. We discuss the results 
of QA employed to reach the ground state of the classical 
spin glass in Sec. IIIII We demonstrate in Sec. IIV1 how 
QA can be utilized to simulate the ground state of a 
quantum spin glass. We conclude with a short summary. 



I. MODEL 

We consider an infinite range Ising spin system whose 
Hamiltonian is given by 



TL = 



N 

E ■ 

;j(>i) 



(1) 



where of is the z— component of Pauli spin, represent- 
ing a classical Ising spin at site i and J^ 's are random 
variables taking up values either +1 or —1 with equal 
probabilities. The above Hamiltonian describes a cluster 
of N Ising spins, each connected to all others through ex- 
change interactions of equal strength (J = 1) but random 
signs. To make the energy extensive in system-size, one 
has to scales the energy with a factor of N 3 ^ 2 , as done in 
the rest of the article. The system is heavily frustrated 
(i.e., no spin configuration can satisfy all the bonds), due 
to the presence of both ferromagnetic and antiferromag- 
netic bonds in random fashion. The high degree of con- 
nectivity (i.e., the infinite range of the interactions) adds 
to the complexity of the problem. For such a system, 
finding the ground state spin configuration for any arbi- 
trary given realization of interactions (the set of Jy 's) , is 
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FIG. 1: Here the complexity of finding the ground state for a typical 30 spin sample of the ± J spin glass (Eq. [TJ is illustrated. 
In (A) we plot the energy versus number of configurations for the sample. In (B) we plot the change in energy of the system as 
one flips the spins one-by-one, starting from a ferromagnetic (all up) state till the mirror symmetric one is reached. In Panel 
(C) we show how energy changes with time (spin flips) for a local downhill optimizer, who flips a spin randomly and accepts 
the move if it reduces the energy. 



known to be an NP-hard problem^. In thermodynamic 
limit the system becomes a non-ergodic spin glass be- 
low some spin glass temperature Tq. To illustrate the 
complexity and the difficulty in finding the ground state 
we consider a typical sample with N = 30 (we stud- 
ied many samples, and the features are qualitatively the 
same as we for this one) as shown in Fig. [TJ The plot 
of the number of spin configurations against the energy 
(panel A of Fig. [T]) of the system clearly shows that there 
is a huge entropy barrier between the ground state and 
the states which has higher energies. A levels with en- 
ergy close to zero has a degeneracy of the order of 10 7 , 
while the ground state degeneracy is of the order of unity. 
This problem could be easily overcome and ground state 
could be reached if the potential energy landscape (PEL) 
had have a more or less monotonic gradient towards the 
ground state, as observed in ferromagnetic samples (or 
may be one having few local non-monotonicity for sam- 
ples with a few disorder and short-range interactions). 
But Here, as illustrated in panel (B), the PEL is quite 
random. In the Fig. QjB), we plot the energy against 
spin flip, as we start from a completely ordered (all up) 
configuration and flip one-by-one all the 30 spins until 
the mirror-symmetric state (all down). In a ferromag- 
net, a similar plot would generate a convex curve (with 
the peak at the zero energy) symmetric about the en- 
ergy axis, reaching the ground state at its both end. But 
here we find a completely random profile far above the 
real ground state energy. There is no helpful gradient to 
guide towards the ground state (if one incorporates moves 
with more than one spin- flips, then the profile becomes 
even more rugged). These two features, that is, the ex- 
tremely low entropy of the ground state and the absence 



of a consistent gradient towards the ground state in the 
PEL makes it absolutely difficult for a random walker to 
find a ground state with the help of a naive local mini- 
mization algorithm, like the principle of going down hill, 
as shown in Fig. [T] panel (C) . Here we show the energy 
of the system evolved using the random single spin-flip 
move, and allowing it to accept the move only if it mini- 
mizes the energy. We find even after 10 7 moves, it fails to 
reach the ground state. One can guess that after getting 
down in energy in first ~ 100 steps, it gets lost among the 
huge number of configuration with higher energies (which 
it does not accept) and cannot find a way to reach the 
sparsely occurring lower energy states. 

The eigenstates of TL (the basis states) are the direct- 
products of the eigenstates of af's. Each basis state 
represents a distinct spin configuration of the system. 
To perform zero-temperature quantum annealing of this 
±J Ising system, we add a transverse field term TL' = 
ft(t) J2iLi °f where erf's are ^-components of Pauli spins 
which introduces probability of tunneling between the 
basis states (classical configurations), and f2(t) is the 
strength of the transverse field. The total Hamiltonian 
is thus given by 

N N 

Htot = H + H'(t) = - J2 ■'■> G ' G , n (*)Z)°f- ( 2 ) 

i,j(>i) i=l 

We start with a high enough value of initially (at 
t = 0) and sample the ground state of TLtot using a zero- 
temperature quantum Monte Carlo algorithm (discussed 
below). During sampling, we reduce the strength Q(t) of 
the transverse field following a linear annealing schedule 



n(t)=n (i-f/T), 



(3) 



where t denotes evolution time. At the end of the simula- 
tion (t = t) we are left with the classical Hamiltonian Ti 
and if r is large enough, the simulated system is finally 
found to be in one of its ground state configurations. For 
low values of r, one generally ends up with a higher en- 
ergy configuration. 



II. THE ZERO TEMPERATURE QUANTUM 
MONTE CARLO METHOD 



we use a zero- 
^ Wede- 



To simulate the ground state of Titot-, 
temperature quantum Monte Carlo technique 
scribe the method to some details here, since it is not 
broadly known, and has been implemented so far only 
to simulate pure systems with short-ranged interactions. 
Here we generalize the implementation for an infinite- 
range system with disorders. 

In this method one makes a linear transformation of 
the form 



W = C1- Ht 



(4) 



FIG. 2: Fhe figure shows the uniform chain of plackets (with 
periodic boundary condition) used to simulate the ground 
state of Titot (Eq. ([2])). A placket (solid circle) in the chain 
is basically the cluster of N Ising spins with a given realiza- 
tion of Jij's, represented by Ti in Eq. ((2j). Fhe interactions 
f([i\, fix+i) between any two nearest-neighbor-plackets are 
determined by the relation VV MaMa+1 = e ~Pf(i*x>i t \+i) ) where 
W is obtained from Titot by the linear transformation (@). 
If the dimension of each placket is d, then the dimension of 
the resulting chain is d + 1. 



be any localized discrete degrees of freedom in general) 
[ii s , as shown in Fig. ©. Each of the /Vs can be in, 
say, p different states. One may note here, that if each 
placket n is a spin cluster embedded in dimension d, then 
the chain is actually a d + 1-dimensional object. Since 
the chain is uniform, its Hamiltonian will be of the form 



where C is a suitable real constant and X is the identity 
operator, such that the matrix representation of W in the 
eigen-basis of Ti is non-negative and irreducible (if such a 
transformation could not be done for an Titot, then this 
method would not be applicable for it). One can then 
consider W to be the transfer-matrix of a uniform chain 
(with Periodic Boundary Condition(PBC)) of classical 
plackets, where each placket is nothing but a classical 
cluster of N mutually interacting Ising spins represented 

by H. 

Now the key point is that one can simulate the chain of 
classical plackets using the elements of its transfer-matrix 
W and in this simulation the equilibrium average of any 
observable (say, energy) related to a single placket is ap- 
proximately equal to the expectation value of the observ- 
able over the dominant eigenstate of W. The dominant 
eigenstate of W in turn, is the ground state of Titot (due 
to the form of the linear transformation between them). 
Thus we actually simulate the ground state properties 
of Titot by simulating the chain. In the next section we 
establish the scheme in details. 



Simulation of a chain of classical plackets using 
Transfer-matrix 



In this subsection we demonstrate that the equilibrium 
averages for a single member of a uniform classical chain 
(with PBC) is approximately equal to the respective av- 
erages (expectation values) over the dominant eigenstate 
of the transfer-matrix of the chain. Let us consider a uni- 
form chain of L identical classical spin clusters (or may 



Ud+i = E f&Xi / iA +i)> 



A=l 

where f(p\, Ma+i) is a p x p matrix whose elements are 
the possible contributions to the Hamiltonian from a pair 
of neighboring spins, as each of them takes up p different 
values independently. The partition function of the chain 
is thus given by 



z = E - E -p 



01 = 1 01 = 1 

P P 



E 



A=l 



g f-jvrn^/ x g-£/(M2,M3) x x g-/3/(/il.01 



E ••■ E ^ 

01=1 01 = 1 



(Ul/i2 X ''0203 



x W u 



where W, 



0X0X+1 



-/3/(0A,0A + l) 



[3 being the temperature 
inverse. Again, since each of fi\ and can take up p 
independent values (i.e., can be hip independent states), 
^0a0a+i defines a p x p matrix W. Hence summing over 



all the indices from /12 to /il and recalling the rule of 
matrix multiplication one gets 



01=1 



Trace{W L 



The matrix W is a transfer-matrix for the chain. If the 
matrix W is symmetric then (it is not the necessary but 
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the sufficient condition) one can write 



where 6 r are the eigenvalues of W ordered by the index 
r, so that \0i\ > \9j\ if i < j. Here a few points are to be 
noted. Since all the elements of W are strictly positive 
at any finite /3, the matrix W is both non-negative and 
primitive (i.e., there exists some finite n, such that W™ 
is strictly positive) . Then according to Perron- Frobenius 
theorem (seei^) the dominant eigenvalue 0\ is strictly 
positive and non-degenerate. Thus 



Z 



E 

r=2 



Here, the leading order error is (02/O 1 ) L and since Oi is 
non-degenerate , 



lim 

L — >oo 







(5) 



for any i =/= 1. 

Now, to see how one can simulate the chain using W, 
one has to note that the probability that the chain be in 
a given state A , in which /ii = iii(A), [i 2 = fJ. 2 (A)... etc, 
is 



P{A) = ( 



x e 



/Z 



= (W Ml(A)M2 (A) x - x W w(A)w(A) ) /Z (6) 

Thus using the conditions of detailed balance, one ob- 
tains transition probability from a state A to another 
state B given by 



P(A -+ B) = 



W 



; M i(^)P2(A) x 



Thus if W is given, we can simulate the equilibrium prop- 
erties (thermal average) of any physical quantity related 
to a placket /i in the chain. To obtain that, we require 
to know the probabilities for the placket /i to be in its 
different possible states when the chain is in equilibrium. 
Let P(/i = k) denotes the probability that the placket is 
found in its k-th state when the chain is at thermal equi- 
librium (at a given f3). If the fc-th state is represented 
by a column vector fc), then these column vectors satisfy 
the matrix relation 

(i\W\j) = Wij, 

where (i\ is the transpose of \i) and the sequence of ma- 
trices implies the proper multiplications between them. 

On the other hand, if \Ei) be the dominant (normal- 
ized) eigenvector of W corresponding to the dominant 



eigenvalue 6\, and if W is hermitian then one can ex- 
pand \Ei) linearly in terms of the basis vectors as 



|£i>=X>fl*>, 



(8) 



fe=i 



where 7^ is the amplitude of the basis state \k) in \Ei). 
Thus in the sampling of \E\) using the basis states |fc)'s, 
the probability of occurrence of the state \k) will be \"f\\ 2 ■ 
Now, one can show that 

P{n = k) = \ 1 1 k \ 2 + 0[{6 2 /e l ) L ]. (9) 

The above equation says that one can sample the dom- 
inant eigenstate \E{) of the matrix W just by sampling 
its basis states (classical configurations of a placket in the 
chain) according to the probability of their occurrence in 
the simulation of the placket at equilibrium in the chain 
(using the elements of W itself, as prescribed in ([7])). 

To prove equation @ , we take any placket in the chain 
and call it n\. Probability that /zi is found in the state 
\k) is 



Pfa = k ) = \ 



M2 H3 ■ 



f*2 V-3 



I <w\ L = W(w) L l*) 

Z [ > kk trace{(W) L } 



(10) 



Above, we have summed up the probabilities of all the 
configurations of the chain, in which /ii = k. Now let 
\6.i) (i = 1,2, ...p) denote the normalized eigenvector of 
W corresponding to the eigenvalue (9^. Then one may 
have a linear transformation between |#i)'s and \k) of the 
form 

i*i>=X>iiwfe> 
A,- 

and the reverse transformation 

ifc} = E^ f )"i^ = E^i^ 

i i 

7 being an unitary matrix. Hence 

w L ifc) = E^fi^> 

=> (k\W L \k) 



E 



using ortho-normality of |#i)'s. Thus, from equation 
pi Ap we get 



P{m = k) 



(k\W L \k) 
trace{W L } 

Si him/Qi) L 

i + E l# i(W) L 
\il\ 2 + o\{e 2 /e 1 ) L ] 
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which proves equation ([9]). 

Thus one can in fact simulate the dominant eigenstate 
of any given suitable (hermitian, non-negative and prim- 
itive) N x N matrix up to a good approximation using 
the above results. One has to define a uniform chain 
(with PBC) of classical plackets, each having TV possible 
configurations. The i-th. state of a placket corresponds 
to the i-th vector of the basis in which the given matrix 
is represented. One then views the given matrix as the 
transfer-matrix for a placket in the chain, and simulate 
the chain using its elements (as prescribed in Q). At 
equilibrium, the probability of getting a placket in its i- 
th state is equal to the modulus square of the weight of 
the i-th basis vector in the representation of the domi- 
nant eigenstate of the given matrix (up to an error of the 
form discussed above). 



(due to the sparsity of the matrix W) by constructing a 
more restricted Markov process to do the samplingii. 

For doing quantum annealing of the same system, we 
start with a high enough value of H and reduce it very 
slowly with time t (Monte Carlo step) following a linear 
schedule. During visiting different plackets in a given 
Monte Carlo step, f2 is however held fixed. The linear 
schedule is specified by Sl(t = 0) = fii„ and The total 
number of Monte Carlo steps executed; Oi„ is linearly 
reduced to zero with t within 95% of the total Monte 
Carlo steps, we thus have f2(f) = Om(l — t/r), where r 
is the annealing time. 



III. RESULTS AND DISCUSSIONS 

B. Implementation of the Monte Carlo 



We now illustrate the implementation of the above 
Monte Carlo scheme by employing it to simulate the 
ground state of Htot given in Eq. ([2]) . Here basis vectors 
fc)'s are the eigenvectors of H, and a classical placket 
is the cluster of N Ising spins with exchange interaction 
described by TL. Now we make a linear transformation of 
the form given in Eq. g]), with C = N(N - l)/2. The 
resulting W matrix is clearly non-negative (since none 
of its diagonal element are all smaller than TV (N — l)/2 
and off-diagonal elements are either or f2(i), which we 
always take to be positive.). Since Htot connects a ba- 
sis state to all other basis states that can be obtained 
by a single spin flip from it, there is no closed subspace 
for Htot- Thus W is also irreducible. It can be shown 
that for a non-negative irreducible matrix, all the results 
of Perron-Frobenius theorem we have used here, holds 
goocU^. Besides, W is of course hermitian. Hence we can 
take W as a transfer-matrix for the chain. It corresponds 
to some interaction f([X\, MA+i) between two neighboring 
(/iA and /iA+i) and some inverse temperature (3 (not ex- 
plicitly important here), given by 

W( AlA , M A+i) = e-' 3/ ^^+ 1 ). 

To simulate the ground state of Htot at a given Q for 
a particular realization of 's, we construct a uniform 
chain of L plackets with PBC. Each placket is a cluster 
of N classical Ising spins (described by cooperative term 
of Htot) connected through the given particular realiza- 
tion of Jjj's (see ©). We start with an arbitrary spin 
configuration (same for all plackets) and a given value of 
VL. In one Monte Carlo step we randomly visit L plack- 
ets. At each such visit we make an allowed move (a 
move whose probability is not trivially zero), such that 
the chain goes from a state A, say, to a new state, say B. 
The probability of acceptance of the move is nothing but 
the transition probability P(A — > B) calculated following 
([7]) (using the elements of W) . While sampling, one can 
easily avoid moves whose probabilities are trivially zero 



We have studied the relaxation behavior of several ran- 
dom Jij samples with N — 30 for linear annealing sched- 
ule (we start with an initial transverse field O m and re- 
duce it linearly with Monte Carlo step, so that it becomes 
zero before last few, 5%, steps. We observe that for an 
annealing of ~ 10 7 Monte Carlo steps, the system reaches 
the true ground state (determined by an extensive search 
method) in almost every case, for a suitably large initial 
transverse field f2j„. We calculate the average exchange 
energy of the chain (over L plackets) in each Monte Carlo 
step, and average that over a few ~ 500 Monte Carlo 
steps. The exchange energy (as given by H of Eq. @ is 
not linear in N and we have to scale it by a factor N 3 ^ 2 to 
obtain the intensive energy density. In thermodynamic 
limit, this intensive energy density approaches the value 
— 0.7633i£ (our finite size results shows some fluctuations 
about that). In Fig. §5§ the relaxation behavior of three 
typical random realizations (Rl, R2 and R3) during their 
annealing are shown. We found that for doing annealing 
of a given sample within a given number of steps, there 
is a suitable range of Sl m . If Sl m falls below the range, 
then the transition probabilities are too low to be able to 
anneal the system within the given time. On the other 
hand, if Sl m is above the range, then the rate of change of 
Q(t) is not slow enough to ensure the convergence to the 
ground state finally (i.e., the evolution in no more adia- 
batic). In Fig. ©, the values of respective fii n 's belong 
to the lower end of the respective ranges. The ranges are 
generally wide enough, and one can find a fl m within the 
range, just by a few trials. 

The relaxation behavior is found to be typically "lin- 
ear" in the sense that the long-time averages decrease 
linearly with time (see lower part of Fig. The re- 

laxation observed in shorter time scale of course shows 
fluctuations around that linear behavior (shown in the 
upper part of Fig. ^ ). This linear nature of relaxation 
is typically seen independent of the details of the partic- 
ular realizations. 
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R2; n in - 7.5 
R3; n in - 2.5 - 
Exact ground-state of R1 - 
Exact ground-state of R2 - 
Exact ground-state of R3 - 
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FIG. 3: In the figure, annealing behavior ol three diflfcrcnt randomly generated realizations (Rl, R2 and R3) ol Jij's are shown for N = 30. 
In Each case the system goes to the exact ground state (shown by respective horizontal lines) at the end of the annealing. In each case 
the annealing time is 10 7 Monte Carlo steps, number of plackets in the chain is L = 600, and each Monte Carlo step consists of visiting 
L plackets randomly and making a random spin-flip trial there. In each case, the transverse field has been reduced to zero from its initial 
value f2i n = 5 following a linear schedule, within the Monte Carlo steps. 
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FIG. 4: In this figure the annealing relaxation (N = 30, L = 600 and Qi n — 7.5) for a particular realization is shown. Here the 
annealing time r = 10 Monte Carlo steps. The upper frame shows the relaxation of intensive energy density with time, when 
averaged over small (~ 500) Monte Carlo steps. The lower frame shows the same relaxation, when the averaging is done over 
a much larger number of (~ 10 4 ) Monte Carlo steps. A linear x 2 ~ n t for the longer time average is shown in the lower part of 
the figure. 



IV. BETTER SIMULATION OF LOW 
KINETIC-ENERGY QUANTUM STATES USING 
QUANTUM ANNEALING 



In a glassy system, where the potential energy land- 
scape has valleys separated by huge energy barriers, sim- 
ulating the ground states (and possibly other low-lying 
states) for low kinetic energy (like the ground state for 
a low value of the transverse of a transverse Ising spin 
glass) using a zero-temperature quantum Monte Carlo 
may be very difficult and time consuming. This is be- 



cause, for small kinetic term, the acceptance probability 
may become very small for higher potential energy states 
energy states, and the system may take a very long time 
to get out of a local potential energy minimum in order 
to visit other equally relevant lower-potential-energy val- 
leys. Thus if one gets stuck is a local minimum far above 
the ground state, at an early, stage of the simulation, 
then it would not be able to reach the low lying valleys, 
whose contributions to the ground state are much more 
significant. This can be remedied to some extent by an- 
nealing the sample quantum mechanically starting with 
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N = 30, Configuration Avreage over 30 samples 
Exact Ground State Energy - 



FIG. 5: Here the variation of the final energy with the annealing time r is shown for the sample size N = 30. Each data 
point is averaged over the same set of 30 disorder configurations. The transverse field Q is reduced linearly from Sli„ — 1.5 and 
the replica number is taken to be L — 600. The horizontal line denotes the exact ground state (averaged over the same set of 
configurations) obtained by an exhaustive search algorithm. 
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FIG. 6: In this figure a comparison between the results of simulation (with and without pre-annealing) of the ground state of the 
Hamiltonian ([2]) for different static values of Q and the corresponding exact results obtained by numerical exact diagonalization 
for the same set of samples. Each data-point represents an average over the same set of 40 randomly generated samples of 
size N — 20. The total number of Monte Carlo steps is 10 5 for each Monte Carlo simulation (including the annealing period 
for the annealed simulations). In the figure, he results of simulations with annealing are seen to be much closer to the exact 
diagonalization results than those with out annealing, for lower values of Q. 



a high value of the kinetic energy, and then reducing it 
slowly up to the low value at which the simulation is 
desired. 

For the Hamiltonian given by Eq. ([2]), simulation of 
the ground state for a small fixed value of the trans- 
verse field n (using the zero-tempcrature transfer-matrix 
Monte Carlo algorithm described here) is found to be 
much closer to the exact result (obtained using exact 
diagonalization^) when the simulation is done following 
an annealing (reducing f2 from a high value to the low 
value at which the simulation is desired) than that done 



directly keeping the value of O fixed to the low value from 
the onset. We compare the results of both kinds of simu- 
lations (with and without annealing), for several random 
samples of the spin glass for N = 20 with the respective 
exact diagonalization results for them (see Fig. [5]). 

We conclude summarizing few points regarding the 
performance of the algorithm described here. The algo- 
rithm discussed here is quite general (applicable to any 
disordered spin system in any dimension) and may be 
used for simulating small system-sizes quite satisfacto- 
rily. However, the moves are very restricted, since each 
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spin flip in a placket requires the two nearest-neighboring 
plackets to be in the same configuration with it. For large 
system sizes, this is too restrictive a condition to move 
freely enough through the configuration space to procure 
a satisfactory sampling rate. In addition, since the accep- 
tance probability for higher potential energy configura- 
tions (like most other zero-temperature quantum Monte 
Carlo algorithms) depends on the magnitude of the ki- 
netic term, it is hard to simulate the ground state for low 



values of the kinetic term. We have shown how quantum 
annealing can be utilized in overcoming this difficulty (at 
least partially). This remedy is is expected to work also 
for other zero-temperature quantum Monte Carlo meth- 
ods in principle. 

We thank M. J. de Oliveira for providing us with the 
basic code for one-dimensional pure Ising chain in trans- 
verse field. A. Das is also grateful to G. Santoro and E. 
Tosatti for useful discussions. 
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